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Effect of Flow on Complex Biological Macromolecules in 
Microfluidic Devices 
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Abstract. Understanding the transport, orientation, and deforma- 
tion of biological macromolecules by flow is important in designing 
microfluidic devices. In this study, epi-fluorescence microscopy was 
used to characterize the behavior of macromolecules in flow in a 
microfluidic device, particularly how the flow affects the conforma- 
tion of the molecules. The microfluidic flow path consists of a large, 
inlet reservoir connected to a long, rectangular channel followed by 
a large downstream reservoir. The flow contains both regions of 
high elongation (along the centerline as the fluid converges from the 
upstream reservoir into the channel) and shear (in the channel near 
the walls). Solutions of X-DNA labeled with a fluorescent probe were 
first characterized Theologically to determine fluid relaxation times, 
then introduced into the microfluidic device. Images of the DNA 
conformation in the device were captured through an epi- 
fluorescent microscope. The conformation of DNA molecules 
under flow showed tremendous heterogeneity, as observed by Chu 
[7,12] and co-workers in pure shear and pure elongational flows. 
Histograms of the distribution of conformations were measured 
along the channel centerline as a function of axial position and 
revealed dramatic stretching of the molecules due to the converging 
flow followed by an eventual return to equilibrium coil size far 
downstream of the channel entry. The importance of shear was 
probed via a series of measurements near the channel centerline 
and near the channel wall. High shear rates near the channel wall 
also resulted in dramatic stretching of the molecules, and may 
result in chain scission of the macromolecules. 

Key Words, microti ui dies, lab-on-a-chip, macromolecular confor- 
mation, hydrodynamics 



Introduction 

Microfabrication processes enable the design and 
manufacture of systems on the microscale. Such 
microsystems are capable of processing and analyzing 
biochemical samples and offer multiple advantages 
compared to conventional systems and protocols. 
Scaling down the system and increasing the number of 
samples processed per chip would result in heightened 
efficiency and reduced cost [1]. The realization of these 
benefits, however, requires understanding fundamental 
physical principles on this scale. Any change in 
conformation or stability accompanying the processing 
of a macromolecule in a microsystem must either parallel 



that which occurs during conventional laboratory testing, 
be benign to the end analysis, or be accounted for in the 
interpretation of data. Toward this goal, this work seeks 
to characterize the behavior of biological macromole- 
cules during flow through a microdevice and to identify 
the critical parameters influencing this behavior. 

A microflow composed of biological macromolecules 
constitutes a multi-faceted problem. Because of high 
surface to volume ratios in a microdevice, viscous forces 
dominate inertial forces and surface forces gain 
significance. Mathematically, Stokes flow and a vanish- 
ing Reynolds number describe the behavior of a 
Newtonian fluid on this scale [2-4]. The addition of 
biological macromolecules to a Newtonian solvent 
increases the complexity of fluid flow by causing non- 
Newtonian effects. For these fluids, viscosity may 
depend upon shear rate, and the behavior of the system 
is a function of the intrinsic properties of the polymer, on 
the solvent characteristics, and on the flow parameters. 
Two quantities predictive of non-Newtonian behavior are 
the polymer relaxation time X, and the Deborah number 
De, a nondimensional parameter equal to the product of X 
and shear rate y: De = Xy. 

Flow through a microdevice is a combination of 
elongational flow and shear flow. Necessarily, the 
predominance of each flow type varies with position 
within the device. For example, shear flow dominates at 
the channel walls, and elongational flow is greatest at the 
channel center near channel contractions. Previous 
studies have experimentally recorded and have modeled 
the conformations of DNA molecules in both a pure 
elongational [5-9] and in a pure shear flow [10-12]. 
These studies reported the stretching of the macro- 
molecule when subjected to either flow as a function of 
De and residence time in the flow. Given these findings, 
we expected to observe similar deformation of DNA 
molecules flowing through a more complex, nonhomo- 
geneous flow in a silicon-fabricated microdevice. Here, 
we focus on the effects of the flow on the conformation of 
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DNA. The effects of the DNA on the flow field are 
assumed negligible; this assumption will be probed in 
detail in a future study. Identifying the conditions 
conducive to deformation and the extent to which 
deformation occurs is relevant to designing a device and 
defining flow parameters that best support the stability 
and functional integrity of the macromolecule during the 
processing and analysis of biochemical samples. 



Materials and Methods 
Device fabrication 

Fabrication begins with a n-type, 4 inch (1 00), single 
surface polished, silicon wafer (see Figure la). Wet 
oxidation of the wafer at 1 100°C for 3 hours produces a 
1 um layer of silicon dioxide on the wafer's surface (lb). 
The first photolithography sequence defines placement of 
the device reservoirs and channel: 8 mm long, 300 jam 
wide (lc-ld). Development of the photoresist then 
allows an oxygen plasma etch of the exposed silicon 
dioxide (le). The remaining silicon dioxide serves as a 
mask for a subsequent reactive ion etch of the device. 
The second photolithography sequence patterns the 




la; n-type Si wafer 




lb: wet oxidation, I u.m S1O2 



Ic: photolithography, mask I 




Id: develop photoresist 




le: plasma etch SiO 2 




If: remove photoresist 



wafer through-holes (1mm diameter) located in the 
device reservoirs (lg-lh). These through-holes function 
as fluidic inlets and outlets and are formed from a deep 
reactive ion etch, or DRIE (li). The application of 10 um 
of photoresist during the second lithography step protects 
the areas of the wafer not etched during DRIE. Removal 
of the photoresist (Ij) permits a reactive ion etch of the 
device features. The depth of these features (reservoirs 
and channel) is a function of etch time, and for this 
device equals 60 (im (lk). Finally, dipping the wafer in 
buffered hydrofluoric acid removes the remaining silicon 
dioxide from the wafer (11). The devices can then be 
separated by dicing with a diamond saw. 

Device bonding and fluidic connection 
The ability of a microscope to resolve DNA molecules 
through a glass coverslip depends upon the thickness of 
the glass coverslip and on the magnification of the 
microscope. An objective magnification of 100 x limits 
the glass thickness to 0.17 mm. In order to effectively 
bond this thickness of glass to a device, a bonding 
protocol using epoxy was developed. First, a spinner was 
used to uniformly apply epoxy to a coverslip. The epoxy 
was then cured at 90°C for 45 minutes. To ash epoxy 




Ig: photolithography, mask2 




lh: develop photoresist 




li: DRIE through-holes 




Ij: remove pliotoresist 




1 k: RIB channels and reservoirs 




11: remove SiO 2 tnHF 



Fig. I. Schematic of the microfabrication process. 
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from the device reservoirs and channel, an oxygen 
plasma was directed through the device via the device 
through-holes. Inspection of the device with a micro- 
scope confirmed the epoxy removal. Polyvinyl tubing 
was then affixed to the through-holes, located at the 
proximal ends of both reservoirs, with epoxy. This last 
step completed the fluidic path, and flow through the 
device could be mechanically controlled with a syringe 
pump. 

DNA visualization 

The biological macromolecule studied in these experi- 
ments was >l-DNA (New England Biolabs). The 
molecular weight of this DNA is 31.5 x 10 6 Daltons, 
which corresponds to 48,502 base pairs. Selection of this 
DNA for experimentation stemmed from two factors: 
First, the contour length of the DNA is long, which 
increases the likelihood for detecting any changes in its 
conformation. Second, within the literature there exists 
substantial documentation of the structure of this 
molecule and of alterations in this structure which 
occur when labeling the DNA with a fluorescent probe 
for purposes of visualization [13-16]. 

As per former researchers, the probe chosen to 
visualize the DNA was YOYO-1 (Molecular Probes) at 
a dye : base pair ratio of 1 : 5. YOYO-1 intercalates within 
the DNA backbone and effectively stiffens and lengthens 
the molecule. Incubation of the DNA and probe occurred 
for a minimum of 2 hours at room temperature. For 
incubation with YOYO-1 and for experimentation, the 
DNA was isoionically diluted in a pH 8 buffer consisting 
of lOmM tris-HG, 2mM EDTA, lOmM NaCl. Due to 
the sensitivity of the dye-DNA complex to light and to 
oxygen radicals, all incubation and experiments occurred 
in the dark and the oxygen scavenger B-mercaptoethanol 
was added to the buffer solution at a final concentration 
of 4% (v/v), respectively. The bulk viscosity of the 
solution was adjusted as indicated by adding sucrose. 

Imaging instrumentation 

For all imaging experiments, an Olympus 1X70 inverted 
epifluorescence microscope equipped with a 100 x / 1.3 
NA Universal Plan Fluorite oil immersion objective and 
a high pressure 100 W, 20 V Mercury burner as the 
illumination source was used. A fluorescence cube from 
Chroma containing the following optical filters was used 
with the probe YOYO-1: excitation 460-500, long pass 
beamsplitter 505, emission 510-560. A high-gain image 
intensifier (KS-1381 Videoscope International, Ltd.) 
coupled the microscope to a high performance mono- 
chrome CCD camera (4900 series Cohu). Data was 
recorded to S-VHS video using a S-VHS Mitsubishi 
VCR. The maximum image capture rate was 30 frames/ 
second. Images were transferred from video to computer 



using a Scion framegrabber and NIH Image software. 
Image analysis was performed on a Macintosh G3. 

DNA rheology 

The rheological characterization of X-DNA was per- 
formed using a Vilastic 3 Viscoelasticity Analyzer from 
Vilastic Scientific, Inc. The instrument produces an 
oscillatory flow in a capillary and measures the pressure 
and volumetric flow rate, allowing the determination of 
both the viscous and elastic response of a fluid sample. 
The Vilastic 3 was interfaced with and controlled by a PC 
containing a Pentium II processor. Data were recorded 
with the Vilastic 3 operating in high drive at a frequency 
of 2 Hz. The drive level was swept from 1 to 10. The 
inner radius and length of the Vilastic 3 measurement 
tube were 0.0508 cm and 6.381cm, respectively. 
Distilled water served as the coupling fluid. A Haake 
K 10 circulator attached to the thermal jacket encasing the 
measurement tube regulated the temperature of the 
sample being measured. All DNA samples were diluted 
as described above and were maintained at a constant 
temperature of 23 °C. 



Experiment 

Characterization of microfluidic device 
The scanning electron micrographs in Figure 2 illustrate 
the microfluidic device geometry in which the DNA 
conformation was imaged. The device shown reflects 
alterations from a previous prototype in that there is a 
single fluidic path (i.e., only one channel) and the inner 
edges of the device are rounded. Experiments demon- 
strate that these changes ease monitoring of flow 
uniformity and facilitate the removal of trapped air 
within the device. For notation purposes, x represents the 
direction of flow, y is the channel depth, and z is 
the channel width. An important design feature is the 
narrowing of the fluidic path from the reservoir into 
the channel which forms a contraction and produces an 
elongational flow. In the present work, we focus on the 
effects of flow on DNA conformation. The effect of the 
DNA on the flow field is assumed negligible and the 
velocity and stress fields are assumed to be identical to 
those of a Newtonian fluid. 

For a Newtonian fluid of viscosity r\ flowing in a 
rectangular channel, velocity and shear stress profiles in 
the y = 0 plane for a channel having cross-sectional 
dimensions of 300 urn x 60 urn appear in Figure 3. These 
represent exact solutions to the equations of motion [17]. 
The solutions for the fully-developed axial velocity 
w(y,z) and volumetric flow rate Q in terms of the axial 
pressure gradient dP/dx are: 
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Fig. 2. SEMS of microfluidic device showing channel entrance and 
reservoir. 



dP\A U^ilf cosh(mz/2^) 1 



to 3 / rfPXf 192<z^ , W2fl 

= 1 ^— > tanh- J — - 

6ti V dxjy n s b fr\ i 5 



I co$(iny/2a) 



(i) 



(2) 



where / is an integer, ~b<z<b and ~a<y<a. To 
obtain the shear stress, the pressure gradient dP/dx as a 
function of Q was determined from equation (2) and 
substituted into the expression for u(y, z). The product of 
the derivative of the velocity u(y = 0, z) with respect to z 
and the viscosity yields shear stress x xz in the y = 0 plane. 
Equations (1) and (2) can be used to calculate the local 
shear rate and De within the microfluidic device. These 
profiles illustrate that the velocity is greatest at the center 
of the channel, and the shear stress is greatest at the wall 
of the channel. In contrast to the parabolic profiles 
characteristic of pipe flow, the high width to height 
aspect ratio of this rectangular geometry results in plug- 
like flow in the central plane; i.e., drastic changes in 
velocity and shear stress occur within close proximity (of 
order 30 urn) to the channel walls. The velocity profile 
along the smaller, transverse direction, u(y, z = 0) is 
nearly parabolic. 



Characterization of relaxation time for X-DNA 
A natural coupling exists between the fluid mechanics 
within a microdevice and the behavior of the biological 
macromolecules flowing through the device. X.-DNA 
exhibits viscoelasticity, and identification of its relaxa- 
tion time X contributes to understanding this coupling. 
Using the Vilastic 3, the viscosity of dilute solutions of 
DNA was measured. This instrument, unlike traditional 
rheometers, employs a piston to force fluid within the 
measurement tube into an oscillatory flow. A nylon 
membrane separates a piston surrounded by silicone oil 
from the coupling fluid (distilled water) in the device 
chamber and from the test fluid in the measurement tube. 
Because the machine lacks moving parts, factional 
effects are avoided. 

During oscillatory flow, the instrument measures the 
pressure gradient along the measurement tube and the 
volume flow through the measurement tube. Knowledge 
of these and additional variables including the operating 
frequency and the measurement tube dimensions allows 
calculation of shear stress, shear rate, and shear strain at 
the tube wall. Determination of the viscosity (or, more 
accurately, the complex viscosity since the flow is 
oscillatory, rather than steady) readily follows from the 
ratio of shear stress to shear rate. 

The relationship between concentration and viscosity 
for X-DNA in a 15 centipoise buffered solution appears 
in Figure 4. The viscosity of the buffer is also shown. By 
fitting the concentration dependence of the viscosity r|, 
the intrinsic viscosity [r\] of the molecule can be 
calculated. The intrinsic viscosity [r|] is related to the 
solvent viscosity r| Jo/v and the solution concentration c 
by: 



n = Tisoiv(i + hk + ^ih] 2 c 2 + ..) 



(3) 



where k { is the Huggins constant. The Huggins constant 
depends on temperature and on the particular macro- 
molecule-solvent system, and typically has a value in the 
range of 0.3-0.5. Defining the specific viscosity as: 



(T1 - Tlsolv) 
^solv 



and substituting for ri^ yields the equation: 
!k =[tl]+A r 1 [ n ] 2 c + -.. 



(4) 



(5) 



The ratio of specific viscosity to concentration is termed 
the reduced viscosity. Plotting this value against 
concentration in the low concentration limit produces a 
y-intercept equivalent to the intrinsic viscosity 
(Figure 5). The nonlinearity of the data in Figure 5 at 
very low concentrations is due to the changing ionic 
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Velocity u(y,z) for 300 urn x 60 \m MicroChannel 
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Fig. 3. Profiles for the fully-developed velocity u{y = 0 ; 2) and shear stress x n in the y = 0 plane in the channel flow of a Newtonian fluid. 



strength of the solution. The ionic strength decreases 
with DNA concentration and the DNA expands in the 
absence of charge-shielding. Taking the linear, higher 
concentration data to represent the behavior of the 
buffered DNA in solution as a flexible coil, the Rouse 
model yields the relaxation time X, as: 



X = 



M]r| solv M 



RT 



(6) 



where M is the molecular weight of X-DNA, R is the gas 
constant, and T is temperature. From this process, the 
intrinsic viscosity was calculated as 3530 ml/g. The 
corresponding relaxation time in a 15 centipoise sucrose 
buffered solvent is then 0.68 second. Since relaxation 
time scales with viscosity, for A,-DNA in a buffered 
solvent at 0.904 centipoise the relaxation time X is 
0.041 second. 



1620 



15.00 




shear me, «*' (k>g«caJe) 



Fig. 4. Viscosity as a function of concentration and shear rate for a buffered solution ofk-DNA. The viscosity of the buffer was adjusted to 
roughly 15 centipoise through the addition of sucrose. 
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Fig. 5. Calculation of the intrinsic viscosity of the DNA in solution. The extrapolation of the high concentration data for the reduced viscosity to 
the y-axis provides an estimate of the intrinsic viscosity [n]. R 2 is the square of the correlation coefficient. 



Characterization of DNA Behavior in a 
Microfluidic Flow 

As a first attempt at understanding how ^-DNA behaves 
in a microfluidic flow, X-DNA molecules were visualized 
at 100 x magnification flowing though the center of a 
microchannel. The molecules were imaged downstream 
of the channel entrance, at a distance approximately one 
third the total channel length from the entrance (Figure 
6). At this location, the flow profile is fully developed 
and lacks entrance effects. The viscosity of the buffered 
solvent was 1 centipoise and the concentration of the 
DNA was 0.30ug/ml. The overlap concentration, c*, at 
which the macromolecules completely fill space without 
overlapping, can be estimated from the rheological data 
as 



Hence the solution under study, at c = 0.001 c*, is dilute. 




Top View of Device fluid out 



Fig. 6. Schematic of device and imaging location for the experiments 
reported in Figures 7-9. 



For this experiment, images of molecules were 
captured for flow rates ranging from 0 ul/hr to 30 ul/hr. 
All other variables were constant. At each flow rate, a 
minimum of 100 molecules was imaged and their 
lengths, as projected onto the direction of flow (*), 
were measured in pixels. Magnification at 100 x of a 
reticle having micron divisions permitted conversion 
from pixels to microns. 

Sample images of X-DNA at each flow rate appear in 
Figure 7. Shown with each flow rate is the maximum 
Deborah number (De) existing within the channel for the 
respective experimental conditions; that is, the Deborah 
number used here is based on the wall shear rate. Since 
imaging of the molecules occurred at the channel 
centerline, where the shear rate tends to zero, the 
observed response of the molecules was not this extreme. 
Because, however, the imaging area spanned about 20 
microns and because focusing with an epifluorescence 
microscope is not exact, the conformation of these 
macromolecules realistically reflects exposure to a low 
shear flow. 

Also shown for reference are fluorescent beads. At 
high flow rates, travel of the molecules may exceed the 
camera capture rate of 30 frames/sec. An optical artifact 
or streaking of the molecules and beads ensues. Imaging 
beads allows detection and quantification of any 
streaking. 

The DNA molecules imaged in Figure 7 seemingly 
elongate along streamlines at increasing flow rates and 
Deborah numbers. Comparison of these images to 
fluorescent beads, however, reveals that streaking 
obscures any stretching as the beads follow the same 
pattern as the DNA. Thus, further attempts to determine 
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Fig, 7. Sample images of DNA and fluorescent spheres flowing through the microchannel as a function of flow rate. The imaging location is 
indicated in Figure 6. The viscosity of the buffer solution is 1 centipoise. 



fluorescent beads 




the effects of flow on the conformation of the 
macromolecules in this solution cannot involve 
increasing flow rates given the existing experimental 
set-up. 

One alternative method for increasing the stress 
imposed on the molecule while circumventing optical 
streaking is to increase the solution viscosity. Repeating 
the experiment in a buffered solvent with a viscosity of 
25 centipoise produces the sample images in Figure 8. 
These images reveal that X-DNA can undergo dramatic 
stretching in a microfluidic flow. Increasing the solvent 
viscosity increased the relaxation time of the molecule, 
and dramatic changes in conformation could be observed 
at flow rates where streaking is essentially negligible. 
Compare, for example, the images in Figure 8 of the 
beads and of the DNA at a flow rate of 20 ul/hour. 

The histogram in Figure 9 contains the percent 
frequency of all recorded molecular lengths. For all 
flow rates, a broad distribution of DNA lengths exists. At 
lOpl/hour, the majority of molecules occupy a coiled 
state. The molecules elongate at 20 ul/hour, and the 
histogram for this data set shifts to the right. By 30 ml/ 
hour, the data histogram looks somewhat bi-modal, and 
the variation in recorded lengths is the most extreme. The 
appearance of a greater number of molecules having 



lengths between 2 um and 3 um at 30 ul/hour relative to 
20 ul/hour may reflect chain scission. 

The molecular deformation observed may be attribu- 
table to a combination of elongational and shear flow. At 
the channel centerline, the shear rates should be extremely 
low and, at the imaging location, the elongational channel 
entrance effects should have partially dissipated. Thus, 
two experiments to elucidate the relative contribution of 
each type of flow to the observed effects logically follow. 
The first experiment examines the elongational flow along 
the channel centerline, and the second explores the shear 
flow at the channel wall. 



Elongational Flow Effects 

An elongational flow exists at the channel centerline; the 
fluid accelerates as it enters the channel from the 
upstream reservoir and decelerates as it exits the channel 
into the downstream reservoir. Therefore, imaging DNA 
molecules moving along the centerline of the device 
lends insight into the impact of a converging and 
diverging flow on molecular conformation and identifies 
the positions within the device most likely to produce 
deformation of macromolecules. By using the fine focus 
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Fig. 8. Sample images of DNA and fluorescent spheres flowing through the microchannel as a function of flow rate. The imaging location is 
indicated in Figure 6. The viscosity of the buffer solution is 25 centipoise. 
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Fig. 9. Histograms of the projected length of the DNA flowing through the microchannel. The imaging location and viscosity are those indicated 
in Figure 8. A minimum of 100 molecules was recorded at each flow rate. 
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/wg. /0. DM4 conformation along the channel centerline as a 
function of axial position in the device. The representative images 
indicate the dramatic stretching of the molecule due to the 
elongational flow at the channel entrance, and the recovery of its 
equilibrium conformation as it travels through the channel and into 
the downstream reservoir. 



of the microscope to locate the top and then bottom of the 
device (approximately 0.6 urn/division), the channel 
centerline was found. As illustrated in Figure 10, the 
imaging position varied along x, the direction of flow, 
while the y and z positions remained constant. The con- 
centration of DNA in buffered solution was 0.3 ug/ml, 
the viscosity was 15 centipoise, and the flow rate was 
5 ul/hour. 

At A, the DNA molecule is in the device reservoir. 
The shape, as represented by the sample images in Figure 
10, is coiled. Having fully entered the channel at B and 
having just undergone a converging flow, the DNA 
molecules are highly elongated. Dividing flow rate by the 
device cross-sectional area gives a mean fluid velocity of 
0.077 mm/sec. within the channel. In this 15 centipoise 
viscosity solution, the relaxation time of the molecule is 
0.682 sec. Therefore, in traveling 0.053 mm along the 
centerline, only about one relaxation time will have 
elapsed, and a large fraction of the elongational stress 
from entering the channel will not yet have relaxed. At 
position C, the approximate midpoint of the 8 mm 
channel, the DNA molecules as seen in the sample 
images are still not entirely relaxed. Greater relaxation of 
the molecule is observed further downstream of the 
channel entrance, at positions D and E. Not until the 
molecule exits the channel at F, however, does the 
molecule recover the conformation observed at A. 

Figure 1 1 is a histogram of the lengths of all DNA 
molecules measured at positions A-F. A minimum of 
100 molecules was recorded for each position. As seen 



with the sample images in Figure 10, the DNA molecules 
elongate in the contraction flow at the channel entrance 
(B). Only after exiting the channel at F do the molecules 
return to a relaxed, coiled conformation observed in the 
device entrance reservoir at A. To ensure optical 
streaking did not bias these measurements, fluorescent 
beads were imaged under identical flow conditions 
experienced by DNA molecules; these data are reported 
in Table 1. 

Since the confidence intervals for the bead lengths 
overlap, streaking of the molecules is not a concern. 

This experiment demonstrates that an elongational 
flow in a microfluidic device can cause DNA molecules 
to stretch, and that the time required for the molecules to 
recover their initial conformation exceeds the relaxation 
time by a large (of order 100) factor. The data suggests 
that an understanding of the fluid relaxation time and the 
detailed kinematics is critical when designing micro- 
fluidic devices. Without this knowledge, the placement of 
too many stresses in series will stretch and degrade the 
molecules and will bias any analysis. 

Shear Flow Effects 

Having examined the elongational flow along the 
channel centerline, the next step is to study the shear 
flow along the channel wall. Figure 12 shows the imaging 
positions for this flow experiment. Molecules were 
imaged at the channel wall and near the channel 
centerline. The ' 'upstream'' wall location immediately 
follows the channel entrance, and the "downstream" 
location precedes the channel exit. At the channel wall 
upstream, the channel wall downstream, and near the 
channel centerline, images were recorded at 5 depths (in 
y): at 6 um, 12 um, 30 nm, 48 urn, and 54 um from the top 
channel surface. The total channel depth was 60 urn The 
top and bottom of the channel were located as previously 
described. In these experiments the DNA concentration 
was 0.3 ug/ml, the solution viscosity was 42 centipoise, 
and the flow rate was 5 ul/hour. Figures 13-15 contain 
histograms of the molecular lengths of DNA projected 
onto the direction of flow. A minimum of 100 molecules 
was measured at each depth. To facilitate a quick 
comparison among depths, the median molecular length 
is reported in the upper right comer. As previously, 
fluorescent beads were imaged to ensure optical 
streaking did not bias the measured lengths of the 
molecules. These results are summarized in Table 2. 

The effects of shear flow on DNA conformation can 
be observed by examining the histograms in Figure 13 
which represent different vertical locations near the 
channel centerline. As indicated in Figure 3, the shear 
rate is highest near the top and bottom walls and 



234 Shrewsbury, Muller, and Liepmann 



1 


1 


position A 




I.I 






<» 9* «« tf W » -.3 « 








position B 




.■ills.. 








a u «# it u u « 






position C 




-III.. 




m ** m *> *s » s» i* *i jf* «» «• «* « *♦ » w 






position D 


.11... 




s* W J* M 


« M V » » Mf «* 


«» SS H » i» W T5 


1 


position E 

li 








position F 









xr molecular length. ujti; y: percent frequency 



Fig. 11. Histograms of projected lengths of DNA on the channel centerline as a function of axial position in the device. The positions are 
indicated in Figure 10. A minimum of 100 molecules was recorded at each position. 
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Table J. Bead lengths at channel center 



Flow Rate 


Bead length in direction of flow (.v) 


0 pl/hour 


1.00 ±0.10 \im 


5 jil/hour 


1.05 ±0.08 \im 



Number of Beads Measured N — 60 



efrannd writ, upstream 



charoifJ wail, downstream 



Direction of Flow 
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channel ccntorfanc ^| 
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Fig. /2. Imaging positions for the flow experiments reported in 
Figures 13-15. 



approaches zero at the channel midplane (y = 30 um). 
This is clearly reflected in increased stretching of the 
molecules near the top and bottom of the channel, and 
relatively little stretching, and a much narrower 
distribution of conformations, near the midplane. Any 
residual deformation of the molecules due to the 
elongational flow associated with the contraction flow 
is modest relative to the deformation induced by shear at 
the walls for these conditions. Theoretically, the shear 
rate at equal distances from the channel midplane is 
identical. Therefore, the distribution of molecular lengths 
at 6 jam and 12 pm should mirror that at 54 pm and 
48 pm, respectively. Differences may arise from dissim- 
ilar surfaces at the channel top (glass) and channel 
bottom (silicon) or from the limited sample size ( > 100 
molecules per depth). 

Figure 14 shows the distributions of molecular lengths 
at different vertical locations near the upstream wall. Due 
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Fig. 13. Histogram ofDNA conformations as a function of depth at the channel centerline location indicated in Figure 12. A minimum of 100 
molecules was recorded at each position. 
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Fig. 14. Histogram ofDNA conformations as a function of depth at the upstream wall location indicated in Figure 12. A minimum of 100 
molecules was recorded at each position. 



the presence of the top and bottom surfaces, the shear 
rate is again highest at these surfaces and decreases — but 
remains high — near the channel midplane. The histo- 
grams at all vertical locations are much more similar in 
shape and breadth than those in Figure 13. In particular, 
comparison of Figures 13 and 14 reveals that the 
molecules are much more stretched, and the histograms 
are much broader, near the midplane of the channel at the 
wall location than close to the centerline. Comparison of 
the two locations at y = 30 um is perhaps the most 
compelling: near the channel center the histogram is very 
narrow and the mean length is 1 .5 urn; at the wall a much 



broader spectrum of lengths is observed and the mean 
length is 1.8 um. 

Figures 14 and 15 reveal a disparity between the 
distributions of DNA lengths at corresponding depths 
near the wall, but at different axial positions within the 
channel. At the downstream wall position, the histo- 
grams are essentially similar at all vertical positions 
within the channel; the hint of a second maximum in the 
histograms in Figure 14 is missing in Figure 15, and 
the population of longer molecules is smaller. In fact, at 
the positions nearest to the top and bottom walls, the 
median DNA length has dropped from 3 and 2.4 um 
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Fig. 15. Histogram ofDNA conformations as a function of depth at the downstream wall location indicated in Figure 12. A minimum of 100 
molecules was recorded at each position. 



Table 2. Bead lengths for shear flow conditions 



Channel position Bead length in direction of flow (x) 



channel wall 6 Jim 


1.01 ±0.11 urn 


channel wall 30 jam 


1.04 ±0.08 Mm 


channel wall 54 u.m 


0.88 ±0.10 jim 


channel center 6 um 


0.99 ±0.09 fim 


channel center 30 jim 


1.05 ±0.09 urn 


channel center 54 jam 


0.89 ±0.08 



q = 5 iil/hr, 

n, = 42 centipoise 



upstream to 1.8 and 1.7 u.m at the downstream location. 
The decrease in median length with axial distance 
suggests that the DNA is undergoing chain scission in 
this high shear flow, and that the molecules are being 
irreversibly degraded. This hypothesis is currently under 
study in our laboratory. 

Conclusion 

The intent of this work was to characterize the behavior 
of biological macromolecules during flow through a 
microdevice. Through microscopy studies of the flow of 
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A.-DNA molecules in a simple microchannel connected to 
inlet and outlet reservoirs, we have observed dramatic 
changes in the molecular conformation due to flow. Both 
the elongational flow near the channel centerline (caused 
by the contraction from the reservoir to the channel) and 
the shearing flow near the walls resulted in profound 
stretching of some of the molecules, and a broadening of 
the distribution of conformations was observed. In 
addition, exposure to extended shearing may lead to 
chain scission — i.e., an irreversible change in the 
molecules. The extent to which the conformation of the 
molecules is affected by the flow is influenced by the 
device geometry, the solution viscosity, the flow rate, and 
the relaxation time of the solution. The Deborah number, 
a dimensionless parameter describing the importance of 
elasticity in the flow, has been used in this initial attempt 
to identify important effects. Higher De flows result in 
increased stretching of the molecules, longer times for 
the molecules stretched along the channel centerline to 
recover their initial conformation, and likely result in 
increased molecular degradation through chain scission. 
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